General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




X-73449-9tl 

i 6 31 11 

TIME-OPTIMAL CONTROL 
OF GRAVITY GRADIENT SATELLITES 
WITH DISTURBANCES 



FRANZ ZACH 


RRvnBER >W9 



GODDARD SPACE FlIGNT COTTER 

6KENKLT. IWtTUfl* 


f 


4, 

Is* ~ ‘ 


/U ^ 


~i 

14 


3f> 


N 

S-, 


V V 


.1 . 



X- 734 -69- 501 



TIME -OPTIMAL CONTROL 
OF GRAVITY GRADIENT SATE LUTES 
WITH DISTURBANCES 


Franz Zach 


November 1969 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 


P?.SCEWNG WW* KAWR NOT «U«t 

CONTENTS 

Page 


1. EOLATIONS OF MOTION FOR A GRAVITY GRADIENT 

SATELLITE (GGS) 2 

3. APPROXIMATION OF GRAVITY GRADIENT SATELLITES BY 

A DUMBBELL SATELLITE (DS) . 5 

4. OPTIMAL CONTROL FOR THE 06 7 

4.1 Optimal Control Laws lor Pitch ( •,) with no 

Disturbance Torques ** 

4.2 Optimal Control Laws for Roll ( ;) and Yaw ( , ) with no 

Disturbance Torques 12 

4.3 Optimal Control Law for Yaw ( .) for I yv = I r/ with no 

Disturbance Torques 12 

4.4 Consideration of Different Reference Values 14 

4.5 Consideration of Disturbances 15 

5. APPLICATION OF THE D6-OPTIMAL- CONTROL TO THE 

GRAVITY GRADIENT SATELLITE 1* 

5.1 Data of the GGS 19 

5.2 Choice of Control Torques 19 

5.3 Comparison of Time Consumed for Damping of librations 

to i - * = . = 0 for GGS's and D6 f s 21 

5.4 Investigations erf the Influence on the Settling Time of 

Initial Conditions not Equal to Zero on all Three Axes ....... 22 

5.5 Pointing to Directions Other Than x * t * 0 22 

5.6 Influence of Disturbances on the Time Consumed for Damping 

of Librations. . . . ; 23 

6. CLOSING REMARKS 23 






PRECEDING PAGE BLANK NOT FfU> 

TIME-OPTIMAL CONTROL OF GRAVITY GRADIENT SATELLITES 

WITH DISTURBANCES* 

ABSTRACT 

The optimal control of spacecraft with gravity gradient coarse stabili- 
zation and active libration damping for precision pointing represents a 
promising approach for advanced experiments, e.g. laser communications, 
where angular position within .001° must be achieved and maintained in 
the presence of disturbances. This paper deals with the time optimal con- 
trol for this case. A subop timai (i.e. nearly optimal) control law for the 
general gravity gradient satellite { GGS ) is obtained by a two-step approxi- 
mation: first the general GGS is approximated by a dumbbell satellite 
(DS) and secondly a small angle approximation is derived. The optimal 
control law for this case is calculated and applied to two representative 
examples of a 24 hour synchronous GGS. The practical results are shown 
and discussed. It can be seen how the optimal control law derived for the 
DS works for the general GGS as a subopt im a! control law. The Maximum 
Principle, which is used in this paper, is also extended to cover operating 
conditions with disturbances. It is shown that disturbances change the 
optimal control law significantly, a result that is very important for practical 
applications of the maximum principle. 

'‘'The work for this paper, was accomplished white the author heldoNafional Research Council 
Postdoctoral Resident Research Associateship supported by the National Aeronautics and 
Space Administration, Goddard Space Flight Center. The author is very much obliged to Mr. 
W fsleyondMr. D. Endres for their valuable suggestions 
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NOMENCLATURE 


The center of all coordinate systems are at the center of mass (c.m.) of the 
spacecraft. All coordinate systems are right handed. 

a ,a .a - gravity gradient coefficient for pitch, roll and yaw respectively 
C - linear viscous damping coefficient of the damper rod 
C 2 - linear spring coefficient of the damper rod 
<i fij - reduced disturbance torque (= T /i ) , i * x, y, z 
DS - Dumbbell Satellite 
GG - Gravity Gradient 
GGS - Gravity Gradient Satellite 

H — inertial angular momentum of the main body of the spacecraft, 
measured in body-fixed coordinates, i - x, y, z 
1 . - moment of inertia of the damper rod about its axis of rotation 

I] - inertia tensor of the main spacecraft body 

r - distance from the center of the earth to the c.m. of the spacecraft 
t - time 

T cM - control torques, i = x, y, z 
T d . — disturbance torques, i = x, y, z 

U - reduced control torques f= j /\ ), i = x, y, z 

1 c i i i 

v - inertial velocity of the spacecraft 

x Q . y 0 . z Q - attitude reference system (x 0 has the direction from the center of 
the earth to the c.m. of the spacecraft, y 0 lies normal to x 0 and in 
the orbital plain such that the component of v in the y 0 direction is 
greater than zero) 

k 

vi 


- coordinate system fixed to the main body of the spacecraft (x^ - 
yaw axis, y b - roil axis, z h - pitch axis) 

- coordinate system fixed to the damper rod 

- pitch angle 

- set value of pitch angle 

- earth gravitational constant 

orbital angular rate, i.e. rotation of the x ri , y n , z 0 system in inertial 
space, measured in the x b , 3^ , z b system, i - x, y, z 

- angles of the damper rod relative to the body - fixed coordinate system 

- time consumed for damping of initial librations 

- roll angle 

- set value of roll angle 

- yaw angle 

> 

- angular inertial velocity of the i*»ain body of the spacecraft, measured 
in the x b , y , : z - system, i - x, y, z 

- orbital angular rate for a circular orbit (= - /r 3 ) 


TIME-OPTIMA L CONTROL OF GRAVITY GRADIENT SATELLITES 

WITH DISTURBANCES* 

1. INTRODUCTION 

The development of control policies for highly precise fine pointing of space- 
craft !s a necessary condition for a number of advanced experiments, e.g. laser 
communications. The gravity gradient technique 1 2 offers a good approach for 
long term coarse stabilization of a spacecraft. However, precise pointing to 
certain angles is required for some experiments, and it may be desirable to 
accomplish this in minimum time. Also for these experiments, disturbances 
such as the solar torque acting on the spacecraft must be counteracted by proper 
control. 

The overall optimal control policy must meet two basic objectives: (1) 
acquisition of the precision pointing mode in minimum time starting from a 
coarse GG mode, and (2) extended term precision pointing by counteraction of 
disturbance torques. This paper will treat only the acquisition requirement. The 
optimal control law as developed herein is applied to two representative examples 
of a synchronous equatorial 24 hour earth spacecraft. The method employed for 
active libration damping is the use of a reaction jet device operating in the micro- 
pound thrust level regime , which can be an ion engine, resistojet., or pulsed plasma 
system. It is envisioned that any use of the ion engine would require special con- 
sideration of control actuations. The main difficulty with the ion engine is that it 
cannot be readily turned on and off for thrusting purposes. One possibility is to 
combine libration damping with orbit corrections. Another possibility is to 


1 


provide a balanced couple which cancels out the orbital effects. The resistojet 3 
has already been studied as a means of optimal large angle maneuvering for a 
synchronous spacecraft. Pulsed plasma systems appear to offer high promise 
for use in both acquisition and holding modes due to their precision impulse bit 
capability and relatively high specific impulse. It is assumed that a minimum 
of six thrust directions is provided or the spacecraft to provide control in pitch, 
roll, and yaw axes. 

The procedure used to develop these control laws is as follows: 

1. The general rotational equations of motion for a GGS including a damping 
mechanism is formulated. 

2. These equations are linearized for a D6 in a circular equatorial orbit 
with small dbrations. 

3. The optimal control laws are developed for this OS. To check the validity 
of these control laws for a more general synchronous GGS, the control 
laws are applied to the simulation of the complete GGS. 

2. EQUATIONS OF MOTION FOR A GRAVITY GRADIENT SATELLITE 

For development of an optimal control law, a GGS which has four rigid GG 
booms and a rigid damper rod is assumed. 

The damper rod will be assumed to rotate in a plane fixed with respect to 
the main body. The spacecraft-damper rod system therefore represents a 
seven degree of freedom system. For pointing the axes of the spacecraft to 
given angles only the equations of rotational motion are of interest for circular 
orbits. These are 
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The torques M ih and M id caused by the gravity gradient effect are. for the main 
body of the spacecraft 
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and for the damper rod 
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P and M sd used in <2) are given by: 




** %<i -( * yd ' ’ 2* 


and 
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»D '“l 1 2 ^2 ' 2 J 


( 4 ) 


where the latter term is the restoring torque exerted on the damper rod by the 
spring-damper mechanism. 


The a t (i = 1, 2, 3) and b (i ~ 1, 2, 3) used in (.3), which give the relations 
of the GG torques to the angular positions, are: 
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The equation of motion for the damper rod is 


t 7 - a yb sin ' l * n .b cos M - *»d 


*,d * *SD 


( 6 ) 


Integration of (1) and (6) furnishes $l %b , U , L ib and 2 . 


The pitch, roll, and yaw angles, which are the Euler angles between the 
body -fixed and the attitude reference system, will be used throughout this report. 
These can be obtained by subtracting l b , the orbital rate which is the rate sit 
which the attitude reference system rotates in the inertial system, from f? b . 
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This gives 


and 





A digital program was written 1 which allows the calculation of , . r and , 

according to Equations (I) through (7>. This program also pays attention to 
elastic collision processes of the damper rod against its stops, which limit the 
rotation of the damper rod in the x — z plane. 


X APPROXIMATION OF GRAVITY GRADIENT SATELLITES BY A 
DOfBBELL SATE LLJTE 

In this section the equations of motion for a OS will be fieri ed from the 
equations developed in the previous section for that of the more general gravity 
gradient satellite (liGS). Thus ft will be assumed that no damper rod exists, 
that I and l rt are the same order of magnitude, and that 



Eq. (1) and (2) mow yield: 


and 


I 

xx 

* xb 

• (1 

yy 

- I ) 

z t 

' y b 

zb 

- T , 

- T „, 


( at) 

I 

yy 

‘ yb 

rl » 

- 1 > 

th 


* 


- V 

(h) 

I 

Z £ 

* 

- ' zb 

1 1 

- 1 ■) 

yy 

" xb 

yh 

- T 

<. * 


• "*> • 

(C) 





9 



Equation (3) becomes: 
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For small angles Equations (9) and (10) can be ap p r o ximated by: 
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where * ~/ r 3 - - 2 and ~ . 

0 ‘ zb -*b =yb 


0 for nearly circular equatorial orbit. 


The sign on the right side of Equation (lib) is changed with respect to Equation 
(9b); this only means that torques afaoot the y axis are considered negative. This 
approach leads to a more comrenien. form for the derivation of the following 
optimal control concepts. 

it is now assumed that 
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as is true, e.g. for the ATS-D and ATS-E spacecraft. Thus Eq. (II) become 


and 
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When d 35 d . * d a * 0 and u , u . , u % are constant each of the Equations (13) 
to (1 5f obncosly gives circles for each of the phase plane plots ), ( t 9 ~ ) or 
( . •:.) respectively. Simulations of the <K38 under the ?sime ccnditxons result 
in circles for the phase plane plots (trajectories) for ttgles ep to 10°. This also 
proves the applicability of the ap p rcKi rngtons used to obtain these eqpatiofft*. 

41 OPTIMAL CONTROL FOR THE DUMBBELL SATELLITE 

Equations (13)* (14), and (13) give a set of decoupled equations for pitch, 
yaw and roll angles of a D,S. These equations make the derivation of an optimal 
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control law using Pontryagin’s Maximum ^nciple practically feasible, ft is 
assumed that the reduced control torques u a , u f and u are bounded in magnitude, 
i.e. 


and 



(16) 


4,1 Optimal Control Laws for Pitch f i > with no Disturbance Torque s 
Let 


and 


Thus (14) becomes 


and 


when d a - 0. 



The theory of optimal control 5 f shows that a vector 



is used to calculate the dot product 
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where 


r*.i 
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' X. 


Next (19) must be maximized by a proper choice of u„ 


Eq. (19) yields 


H P, * 2 p 2 
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which is maximized by 


U. K„ sgn p. 


The theory of optimal control shows further, that 
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Solving (21) gives 
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where A, r - 0 are integration-constants Eq. (29) can now be written as 



The solution of (18) is then: 


and 


where 



( 24 ) 
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x 20 - x 2 ( t - 0) - ( t - 0) 
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Stcan be seen that x, ( t) and x 2 (t)/ |TaJ form a circle in the x t (t) - % 2 ft)/ - 


plane. 

0 

The center of this circle lies at 

K 

x — for u r K 
a a *■ 

a 

and at 

K a 

* = - - for u. = . 

<*■ 

To find the optimum switching law we first consider initial conditions (x J0 , 
x 20 ), which lie on a circle which passes through the origin. Comparing Eq. (23) 
and (24) , it is clear that optimal trajectories cannot stay longer than for an 


angle yT t - ■ on a specific circle. Therefore one can obtain the geometric 
locus of all points with initial conditions (x , x /fT - 0), which can be transferred 
immediately into the origin, to be the two semi-circles shown in Fig. 2. 

Now assume initial conditions which require at least two switching points 
according to Eq. (23). Also assume the last trajectory will lie on 3_, . Eq. (23) 
and (24) show , that the next to the last state with u - * K , must last for 

a 



It can be seen further, that this state is represented by a semicircle with its 
origin at K a / a i , with its endpoint on S_j and its starting point at a semicircle 
with its origin at 3 K a /a and a radius of K a /a a (see Fig. 3 and Ref. 5 and 6). 
When we apply this procedure for all such states leading to S_ t , we see that the 
switching circle for switching from - K a to + K a is given by another semicircle 
with its center at 3 YiJ \ and a radius of K a /a a . When we continue this procedure 
for all initial conditions in the x x - x 2 /y§“ plane, we see that the switching 
circles are located as in Fig. 3. 

Therefore the centers of the switching circles are located at 

K 

' ( 2n - 1 ) -2- . n - • • • - 2. - 1 ,.0. 1 . 2 - - . £5) 

a . 

where n is the number of the switching circle. 
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4.2 Optimal Control Laws for Roll ( ;) and Yaw ( with no Disturbance 
Torques 

When Eq. (13), (14), and (15) are compared, no difference in their form is 
noted as long as a ¥ 0, a „ ¥ 0 and a ¥ 0. This is valid for a spacecraft with 
different principal moments of inertia. Therefore for this case the same optimal 
control law can be applied to Roll and Yaw as is derived for Pitch in Section 4.1. 
This is obtained by substituting (a a . K a ) by (a A ,K. ) or (a . K ) respectively. 
However, another control law must be derived, when two or three of the principle 
moments of inertia are equal. This will be done in the next section for yaw, when 

1 = 1 . 

y y * z 

The flow chart of a digital program for the control laws developed in this 
section and section 4.1 is given in Fig. 4. In this program a deadband is defined 
as 



with a = a ,a . or a respectively. 

i a - f 

Control is switched off completely, when the state of the system is inside 
this deadband (loop 1). Disturbances or later the application of the exact GGS — 
program can drive the system ot& of the deadband again (loop 2). 

4.3 Optimal Control Law for Yaw (f ) for L = I with no Disturbance Torques 

For a symmetrical spacecraft I yy = 1 21 » I nx and hence a b becomes 0. For 
this case another optimal control law must be derived for yaw {$): 
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Let 


and 


then (12) becomes 


and 


1 * l 
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X 1 X 2 
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Using (16) and (26) the Hamiltonian becomes H s x 2 + p u , , which is maximized 
by 


u = 


X sgnp 2 . 


Applying 


dp, _ .|| 

— • — (' i - 12) 

dt r)X. 


again gives 
and 
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which results in 
and 




P, -0 

p 2 = 


p, =B 


'I 


P 2 = - Bt 


* C • .j 


From Eq. (27) and (29) one obtains 


(28) 


(29) 


and 


t 2 

Xf ~ - sgn (C - Bt) + Dt + E 


x 2 t sgn (G - Bt> * D 


(30) 
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Thus all trajectories are parabolas. From (29) it is evident, that p 2 can change 
sign at most once. Therefore the complete optimal trajectory from a given initial 
condition in the x 1 - x - plane consists of one or two segments of a parabola, 
where the second parabola leads to the origin. It is evident, that the parabolas 
leading to the origin are also the switching lines (see Fig. 5 and Ref. 5 and 6). 

For initial conditions with x 20 = 0 the time consumed for reaching the origin 
can be calculated to be 


- 2 

Figure 6 shows the flow chart for a digital program of this control law. 

4.4 Consideration of Different Reference Values 

The optimal control laws in the above sections lead a , and --P to the origin. 

We will now consider the case for reference values (set values) a a , -p ,v>- 4 0. 

* 

Considering Eq. (13) through (15) with a , a. , a. =/ 0 one sees that the follow-. 
mg new equations can be written, which contain an additional control term of thrust 
to balance the gravity gradient torque at the set value: 



and 


• • 

* + a -i lj V + a .A ~P % • 

± K^ + a* ^ 

'i r \ a- - X K a + a a • 


But Eq. (31) also can be written as 



• • 

P i + % Vj = f % * 

•• , 

" • ± 
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and 
where 

and 

a, = a - a . 

1 S • 

Therefore Eq. (17) - (24) can be applied as before, now leading and a, to 0 
and therefore > , > and a to the desired setpoints. In Fig. 3 the switching lines 
will shift and also the optimal trajectories. This results in Fig. 7. 

A similar approach is valid for yaw (y) when a , - 0. The switching lines 
of Fig. 5 must simply be shifted so that they lead to y s . 

4.5 Consideration of Disturbances 

The spacecraft usually is subjected to various disturbances, which are 
different in magnitude and time history. However, there are two typical kinds 
of disturbances, which should* be considered: a constant bias torque and a dis- 
turbance torque with a 24 hour period. A constant disturbance torque can be 
caused e.g. , by attitude thruster misalignments: station keeping thrusters can 
introduce disturbance torques where the thrust vector is not aligned through the 
center of mass. A disturbance torque with a 24 hr. period is caused, e.g. by 
solar pressure acting on the synchronous spacecraft. 

First a constant disturbance torque will be considered. 

a + a . a - u + d ’ , 

a. ad 


Eq. (IS) now is 

( 32 ) 


*i ' \ * K « ■ 


r l 
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where d a is the constant disturbance torque divided by I zz . There are two 
approaches to solve this problem theoretically. First d a could be measured by 
observing the shift of the trajectories and then counteracted by a control torque 
of the same magnitude. The total reduced control torque would be now: 


u 


C Z 


z z 


K a - <*. • 


(33) 


This would lead to the same switching lines for optimal control as in Fig. 2, 3 
and 7. However the practical hardware constraints show that control torques 
according to Eq. (33) are much more difficult to implement than control torques 
of the form 


Therefore a modification of the switching lines of Fig. 2, 3 and 7 is found as 
follows: 


From (32) one proceeds as was done for (17) - (22): 

X| - x 2 , 

k - a x + u + ci 

2 a 1 a a 

and 

* 

H = p r x 2 + p 2 ( - a a Xj + u o + d a ) 
which again results in condition (23). 
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Substituting ± by ± K a + in Equation (24) shows that ail trajectories are 
shifted by d x 7 a, in the direction of the x -axis and that the radii change. Similar 
considerations which led to Figure 2 and 3 give Fig. $ for the optimal switching 
lines. 

Similar considerations hold for optimal control of .. when a =0. We only 
have to substitute ± K by ± K. +- d in Eq. *30). The equations for the switching 
lines read now: 



ft is obvious, that the above calculations are only valid lor R, • d for i= ; 

For slowly changing disturbances optimal control could be implemented by 
changing the switching lines according to Fig. & and according to the momentary 
value of d t . The interactions of changing disturbances and c h a n ging of switching 
lines bring up serious problems of identification of disturbances and problems 
of stability, which are presently under investigation. Two approaches to optimal 
control while the disturbance torque is varying are shown in Fig. 9. It is assumed 
that the optimal trajectory starts without disturbance at x , » 20 /y» and reaches 
S* and P t and follows S* {after switching control from -Kto +Kj until a dis- 
turbance d occurs at the instant the state is at P 4 . The trajectory now changes 
from that of S^. to , where ^ depends upon the disturbance. The S, + tra- 
jectories are circles with their origin shifted by d /& is direction of x . Two 
procedures bow exist: 

* S , rvploee* Sg, S_ f«p*eees 5_-j.fR smctim 4 1 
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a. Retention of ♦ K At P^ , thereby following trajectory R u to P 3 and there 
switching to - K. When S f „ m reached at V 4 , control is switched bock to + K and 
S l+ followed to the origin, or 

b. After occurence of the disturbance d at P , control Is switched immediately 
to -K, until S dit is reached at P, » S r}+ then is followed to the origin by changing 
the control bock to * K. 

However, for the further investigations of the optimal control law in this 
paper it is assumed, that T d ' K for all three axes. Thus, it is not necessary 
to apply the refined control taw given in Fig. 8 and 1 The approach to the origin 
in presence of small disturbances can be seen from Fig. 10. 

ft can be seen from Fig. 10, that the positive x axis may serve as an ad- 
ditional switching line when the trajectory intersects from the negative % 2 direc- 
tion. The negative x { - axis forms the additional switching line in die case when 
the trajectory intersects coming from positive % -values. A more sophisticated 
control law in uie presence of disturbances is suggested for future work in this 
area as follows: switching from + K to -K or vice versa shall be implemented, 
when the distance from the set point (here origin) starts to increase. This 
should at least give stable performance under all working conditions. 

5. APPLICATION OF THE DS-OPTHIAL-CCWTROL TO THE GRAVITY 
GRADIENT SATELLITE 

This section will discuss the results of app ly ing the D6 control to the GGS 
computer simulation. Two GGS - co nfi g u rations are used. The first, which will 
be called GGS 1, is characterized by 1 = I > I and no damping mechanism, 

yy tut ** 

the second (GGS 2) is characterized by the moments of inertia of the ATS -D satellite. 


5.1 Date of the GGS 

a. For the symmetric case (GGS 1) 


1000 0 


0 


0 48000 0 


0 


0 48000 


( tn si ft 2 > 


I . 0 and C f C. 0 

f 12 


b. for the uo»yxnrnetric case (GGS 2) 
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5.2 Choice of Control Torgoes 

hi order to choose the proper magnitude of control torqpe it shooid be noted, 
that a necessary condition for controllability of a gravity gradient s at ell ite is, 
that in the whole pointing region for each axis 


T ^ T 
■e «g 




m 
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where T is the control torque exerted by the ion thruster, T is the restoring 
torque due to the gravity gradient effect for a particular set of value and is 
the disturbance torque. To point in a particular direction . . ; . 1 a certain T 
has to be counteracted; this torque can be calculator by Eq. (10). Since T 
changes only when the set values of •• . ; , and . are changed and therefore mi»ch 
'ess often than T , , counteraction of T by a control torque of the same magnitude 
can be assumed. Pointing is desired to every point on the earth disc as seen from 
a synchronous spacecraft. Therefore -8.5" _ , ; . 8.5" is valid, and a maximum 

gravity gradient torque in this region for the pitch ( • > - axis of the asymmetric 
spacecraft is given by 


A disturbance torque i assumed to be caused by solar pressure with a peak value 
of 10 1 ft *. Eq. (35) now yields 

T - 35 * KT* ft lb 


To give the same T /I ratios on the other axes of the asymmetric spacecraft 
and on all axes of the symmetric spacecraft, the following torques are chosen 
for the asymmetric spacecraft 


and 


T. t r 28. 5 * 10~* ft lb 
T *±66* l O' 6 ft lb 

C**.' 


and for the symmetric spacecraft; 


T - ? 10" 4 ft lb . 

Ctt 

T ci = 1 l V* ft lb 
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and 


T„ - 2 09 • 10 * ft lb 

5.3 C omparison of Time Consumed for Damping Librations to 
■ - ; - , =• 0" for GG Satellites and PS 

The optimal switching lines given in Fig. 3, 5 and 7 are applied to all three 
axes simultaneously. To compare the results obtained for a GGS with the ideal 
case of a DS, the time consumed for reaching an angle of 0 - after a certain initial 
condition is plotted in Fig. 11 oik! 12. For Fig. 11 the time consumed for a DS is 
calculated by' applying the optimal control law derived in Section 4 to the equation 
x a ax = u where a - a, and u = ± K., . Two different , called and a a<t are 
used for GGS 1 and GGS 2, respectively. Therefore two different DS must be used 
ior comparison. DS 1 uses a , whereas DS 2 uses a . . The lower time consumed 
by a GGS for higher initial conditions can be explained as follows: the coupling of 
all three axes transfers some of the energy of the pitch axis to roll and yaw axes 
where optimal control is also being applied. Therefore energy is taken out of 
the whole system faster than for the uncoupled case of a DS. The region of 0 be- 
tv 15" and 60" for GGS 2 shows on the other hand, that this effect is outweighed 
by the less exact approximation of an asymmetric GGS by a DS. Therefore a 
slightly higher time is consumed for this region than is predicted. 

Figure 12 shows the equivalent comparisons for initial conditions 4 0 only on 
roll (; 0 ). Approximately the same a ± is used for both GGS 1 and GGS 2. Therefore 
approximately the same equivalent DS is valid for comparison of both GGS. 

Similar explanations for the shape of the plot hold as for pitch. 
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Since the yaw angle is not significant for pointing to certain points on the 
surface of the earth, plotting of the graphs for yg v equivalent to Fig. 11 and 12 
are not given here. It is just mentioned, however, that for GGS 1 consistency 
with predicted values is achieved, whereas for GGS 2 some deviations yf predicted 
values are observed. An exact explanation of this problem would involve a dis- 
cussion of the exact equations (l) to (7) and goes beyond the scope of this paper. 

5.4 Investigations of the Influence on the Settling Time - of Initial 
Conditions not Equal to Zero on ail Three Axes 

To show the influence on - of other initial conditions than in Section 5.3 
; 0 = , n ^ 0 and i o = o 0 - 0 are chosen. Fig. 13 shows - for pitch 
for GGS 1 and GGS 2 in comparison to the OS. Fig. 15 presents the graphs for 
roll. For yaw the same remarks are valid as in Section 5.3. 

5.5 Pointing to Directions Other Than a = ^ = 0° 

In the previous investigations, it always was assumed the pointing should lead 
to.i s = = 0 n . Now pointing to \ with -8.5 9 _ i n , < 8.5° is considered. 

The simulations show that application of the control law derived in Section 4.4 for 
reference values 4 0 gives no difference in r when pointed to set values a a , 
other than a a - \ - O'. This assumes that \ or - + % now has the 
same value as - 0 or x Q respectively in Sections 5.3 and 5.4. This assumes further 
that the same control torque T is still available and the torque necessary to 

c 

counteract the gravity gradient torque T ^ at the set value is supplied in addition 
to T r . ~ is increased if counteracting of T, r is provided by control thrust while 
decreasing the thrust T. available for actual control purposes. Fig. 15 gives ” 
for the OS as a function of K = f c /I. The parameter 0 is the difference between 
initial Condition and the set vt.*ue. 
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5.6 Influence of Disturbances on the Time Consumed > for 
Damping of Librations 

Investigations concerning control under the influence of a sinusoidal disturbance 
torque with a 24 hour period were carried out. The peak value of T (J is assumed 
to be 10 5 ft lb. * is plotted in Fig. 16 for pitch ( i 0 4 0, . 0 * . 0 * 0) and in 
Fig, 17 for roll ( - 4 0, .■ = , A =0) with V - l = 0 in both cases. 

Two phases for the disturbance torque are chosen, the first is 0° and the 
second is 180° relative to the starting point of control application. 

It can be seen that under favorable phases of the disturbance torque the effect 
on control action nets in no increase of ' . One can also ccnsider cases, where 
- could be decreased; whereas it can be seen from Fig. 16 and 17, that under 
unfavorable circumstances * is increased by 50%. 

6. CLOSING REMARKS 

It has been shown, that the time- optimal control law for a dumbbell satellite 
can be applied to a gravity gradient satellite as suboptimal control. 

The next major consideration must be development of a control policy for 
the extended term precision pointing mode. Of special interest is the formulation 
of an adaptive disturbance torque computer for operation in a scallop mode limit 
cycle. It is also intended to expand the performance index for the acquisition 
mode to include both time and fuel weighting factors. Methods of state determina- 
tion will also be assessed during follow-on application studies. 


23 




REFERENCES 


1. Barrett, C. C., "The Development of a Mathematical Model and a Study of 
One Method of Orbit Adjust and Station Keeping Gravity -Oriented Satellites," 
TND-3652, Nov., 1966, NASA. 

2. Roberson, R. E. , "Attitude Control of Satellites and Space Vehicles," in: 
Advances in Space Sciences, Vol. 2, New York: Academic Press, Inc., 1960, 
pp. 351-436. 

3. Isley, W. C., "Optimal Control Application for Electrothermal Multijet Sys- 
tems on Synchronous Earth Spacecraft," Journal of Spacecraft and Rockets, 
Vol. 5, No. 12, Dec., 1968, pp. 1444-1451. 

4. Frick, R. H., "Perturbations of a Gravity Gradient Stabilization System/’ 
Report AD 623279, The Rand Corporation, Santa Monica, California, Sept., 
1965. 

5. McCausland, I., "Introduction to Optimal Control," New York: John Wiley 
& Sons, Inc., 1969. 

6. Jones, G. S., Strauss, A., Zvego, G. P., "Mathematical Aspects of Control 
Theory and the Problem of Satellite Attitude Control," NASA Contract NAS 
5-9172, Final Report, January, 1966. 


24 



gg— GRAVITY GRADIENT BOOM 


Figure 1- Gravity -gradient satellite and y b lie in the plane of the four booms and 
together with establish a right-handed coordinate system fixed to the main body of 
the spacecraft x d , y d , r 4 establish a right-handed coordinate system fixed to the 
damper rod. 





Figure 4— Simplified flowchart for trajectories with optimal control according to Fig. 3 for 

gravity gradient coefficients a a , a^, a ^ * 0; 
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Figure 7- Sw itching lines aM optimal fra feet or re* far *j - a, £ cr * and *2 r ^ ~ ^ ^ ***!*•€ 
fively. Also fa, 1C) = fa,, KJ. (a^, K^) or -{* b , K j) re*peet ively 





Figure 9-Optimal trajectories under the influence of changing di sturfaonces The opfimot 
trojectory starts at Sjg >a) without disturbance. A disturbance d uccurs when the 

state is at P 2 The switching lines 5* and S_ are then replaced by and 5^.. respec- 
tively. Two possible control procedures that lead the state to the origin we: (a) a policy 
that yields the path Pj, Pj, P 4 , origin, or fb) a policy that yields the path Pj, P 3 , P 4 , 
origin. For K, a, d, *j and sct P ig- 8 - 
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(«q/ deg 


Figure 11 -Comparison of time consumed (r) for damping of initial pitch angles (a^) to 0°. 
The initial conditions on roll (£g) and yaw (0 q) are 0°, x^ - ~ "0 

"V GGS 1 " GGS with two principal moments of inertia equal 
♦ GGS 2 - GGS with all principal moments of inertia different 
OS 1 - Dumbbell satellite with a^and K a as GGS 1 
DS 2 - Dumbbell satellite with a a and K a as GGS 2 





Figure 12 -Comparison of t for damping of initial roll angles (<^)to0°. The other initial 
conditions are: Og 0, Xg g - \£g ~0. o^ and K ^ ore equal for both GGS and there- 

fore the same D5 is applied for comparison. 
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i»q/DEG 


Figure 13-Comparison of r for damping of initial pitch angles (oq) to 0° considering initial 
conditions on the other axes not equal to 0°. the other initial conditions are = ■Pq = 

DS 1 again is for comparison with GGS 1, DS 2 for comparison with GGS 2. 



0 


4 8 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 

(fy/DEG 

Figure 14-Comporison of r for dampingof initial roll angles (j^) to 0° with initial conditions 
on the other axes not equal to 0°. These other initial conditions are Og = Og = Xg = 

f Q - 0. Both GGS'1 and GGS 2 can be compared with the same DS. 
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Figure 15-Time consumed by o DS for changing the direction of pointing by 0 

8° as o function of the reduced torque K =T 'I. 


T , 6°, and 






oq/ DEG 

Figure 16-Influence of disturbances on r for initial pitch angles for GGS 2. 
i'Q ; %j ~ ^0 ^0 The disturbances are sinusoidal with a peak 

value of 10 -5 ft lb, a period of 24 hours and a phasing of 0° relative to the 
start of the control in the first case and 180° in the second case, 
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Fr§*r* W-\tA\wmc» «f ^ifwWwcw an r for rnttwl rail M|Im far CCS 2. 
vp-'&t =0/ . 2 = ®' Tb AuwWnc** ara <b wa a» far FSf. 16. 




